Spin relaxation of two-dimensional holes in strained asymmetric 

SiGe quantum wells 

B. A. Glavin 

V. E. Lashkar'ov Institute of Semiconductor Physics, 
Pr. Nauki 4I, Kiev 03028, Ukraine 

K. W. Kim 

Department of Electrical and Computer Engineering, 
North Carolina State University, Raleigh, NC 27695-7911, USA 

Abstract 

We analyze spin splitting of the two-dimensional hole spectrum in strained asymmetric SiGe 
quantum wells (QWs). Based on the Luttinger Hamiltonian, we obtain expressions for the spin- 
splitting parameters up to the third order in the in-plane hole wavevector. The biaxial strain of 
SiGe QWs is found to be a key parameter that controls spin splitting. Application to SiGe field- 
effect transistor structures indicates that typical spin splitting at room temperature varies from 
a few tenth of meV in the case of Si QW channels to several meV for the Ge counterparts, and 
can be modified efficiently by gate-controlled variation of the perpendicular confining electric field. 
The analysis also shows that for sufficiently asymmetric QWs, spin relaxation is due mainly to the 
spin-splitting related D'yakonov-Perel' mechanism. In strained Si QWs, our estimation shows that 
the hole spin relaxation time can be on the order of a hundred picoseconds at room temperature, 
suggesting that such structures are suitable for p-type spin transistor applications as well. 

PACS numbers: 73.21.Fg, 71.70.Ej, 71.70.Fk 
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I. INTRODUCTION 



Recently, there has been considerable interest in the field of spintronic applications. They 
range from various approaches in logic and memory elements to quantum computation.— 
Many spintronic device concepts rely on the so-called Das-Datta spin transistor,— where 
the transfer of spins between the spin-polarized source and drain contacts is controlled 
by the gate bias through the perpendicular (i.e., confining) electric field in the quantum 
well (QW) channel. Experimentally, the efforts have concentrated mainly on n-type III-V 
devices; in particular, those with InGaAs QWs. The reasons for this preference are (1) 
easy control of electron spin states in this material^ and (2) availability of III-V diluted 
magnetic semiconductors.— Much less attention has been paid to the hole-based spintronic 
applications, whose rationale actually originates from the considerations related to the bulk 
semiconductors. Indeed, in bulk cubic semiconductors the top of the valence band at the 
Brillouin zone center is four-fold degenerate (including spin), and hybridization of the hole 
states for finite wavevectors depends on the direction of the wavevector. As a result, hole 
scattering is supplemented by simultaneous change of spin, and the so-called Elliott- Yafet 
mechanism of spin relaxation is very effective.— In fact, the hole spin relaxation rate is close 
to the momentum relaxation rate; i.e., nonquilibrium spin relaxes too fast to allow control 
in any realistic device. This picture is not exactly correct in the case of two-dimensional 
(2D) holes confined in a QW, where the degeneracy at the zone center is partially lifted 
due to quantization. Moreover, the hole spectrum in an asymmetric QW is completely 
nondegenerate for nonzero wavevectors, qualitatively resembling that of electrons. The 
effective Hamiltonian for 2D holes in GaAs QWs were analyzed, for example, in Ref. 0. 

In this paper, we address the problem of the hole spectrum and spin relaxation in strained 
SiGe QWs. Previously, spin splitting in such QWs was estimated based on the observation 
of the circular photogalvanic effect^ and transport measurements.— In Ref. |7|, possible sources 
of spin splitting proportional to the in-plane hole wavevector were examined briefly. Here, 
we present quantitative calculations of spin splitting based on the Luttinger Hamiltonian 
approach. We take into account the biaxial strain inherent in the SiGe heterostructures 
which gives rise to stronger splitting of the hole subbands. We consider both the cases of 
tensile and compressive strain resulting in the light-hole (LH)-like and heavy-hole (HH)- 
like ground state, respectively. Based on the results for spin splitting, we analyze the 
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spin-relaxation process for 2D holes. For asymmetric QWs, it is found that the D'yakonov- 
Perel' mechanism is more important than the Elliott- Yafet mechanism. According to our 
calculations, the hole spin mean-free path in a strained Si QW can be as large as a micron at 
room temperature. This estimate suggests the feasibility of p-type spintronic applications 
based on Si, particularly in view of recent advances in group-IV magnetic semiconductors.— 

II. BASIC EQUATIONS 

First, we need clarify what we mean by the term "spin". Due to the degeneracy at the 
zone center, the spin-orbit interaction leads to strong hybridization of the hole states and 
the total angular momentum must be considered. If some asymmetry is introduced such as 
the confinement in a QW, the degeneracy between the LH and HH states is lifted at least 
partly. When the resulting doublet is well separated energetically from the other states, it 
can be treated by using a spin Hamiltonian with the effective spin 1/2. So, in speaking of 
spin states, we mean the states of this quasi-degenerate doublet. 

Let us describe the model used in our calculations. We consider the case of a QW grown 
along the [001] direction that is subject to biaxial strain. We start from the following 6x6 
effective-mass Hamiltonian: 

if = i7f +i7P + if, + f/(z)/6. (1) 

Here = + i^^^"^ is the Luttinger Hamiltonian with H^^^ corresponding to the part 
with kx^y = 0, is the contribution due to the biaxial strain, U{z) is the confining potential 
which forms the QW, and Jg is the 6x6 unity matrix. The explicit expression^ for H is 
not provided since it is readily available in the literature (see, for example, Ref. lid ). Note, 
however, that in the following the top of the HH band is used as the reference energy. As 
in k ■ p theory, if}^"^ is treated as a perturbation. Following the conventional notation, we 
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choose the zeroth-order (i.e., unperturbed) wavefunctions as 
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where the plane- wave factors exp[i(A;a;X + /cy?/)] are omitted for simphcity. In Eq. the first 
subscript of \l/ (i.e., 1 or 2) denotes the doubly degenerate states for kr^.y = 0, n is the subband 
number, and /, /i, or s represents the LH, HH, or spin-split (SS) states, respectively. In fact, 
if the strain is strong enough such that the corresponding strain energy is comparable to 
the spin-orbital gap A in the bulk material, the so-called LH and SS states are hybridized 
even for k^^y = 0. Hereafter, we term the states "LH" and "SS" which take the genuinely 
light and spin-split nature in the limit A — > oo. This hybridization is evident from the 
eigenvalue equations for the LH and SS envelope functions x- 
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Here x stands for I or s, A and B are parameters of the hole spectrum,— ttiq is the free electron 
mass, I2 is the 2x2 unitary matrix, and the deformation energy mentioned previously 
is = h{2uzz — Uxx — Uyy), where b is the deformation potential constant and u is the 
strain tensor. The envelope functions for HH states xi!^^ determined by the conventional 
Schrodinger equation with an effective mass m^'^^ = mo/(v4 — B). 

In general, further steps require numerical solution of Eq. 1^ to obtain the spectrum and 
envelope functions of LH and SS states. Considerable simplification is possible, however, if 
the quantization energy is much less than or A, which is often the case for strained SiGe 
QWs. Under this approach, it is adequate to use an "effective mass" approximation for the 
solution of Eq. Q, which treats the nondiagonal part of H2 containing the 2;-derivatives as 
a perturbation. Subsequently, we obtain 
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Here the spectrum -E^''**-* and the envelope functions Xnu Xns are determined by the conven- 
tional Schrodinger equation with the effective masses 
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In Eq. (jH), tij are the elements of the unitary transformation matrix T that approximately 
diagonalizes H2 (for U = and d/dz = 0): 
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where E^'"''^^ are the LH and SS states in a strained material with no confinement (i.e., U = 0): 
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These basis functions result in nondiagonal terms in the Hamiltonian H. Using the 
perturbation method,-^ it is possible to transform it to the quasi- diagonal form in any desired 
order in H^^\ Under this approach, the effective Hamiltonian for the quasi-degenerate 
ground subband can be written as 
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where a are Pauli matrices. Restricting ft to the terms proportional to k and and 
neglecting the corrections to the longitudinal effective mass, we obtain for the ground LH- 
like state 
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In Eqs. ffTTj) and (fT^. we use the following overlap integrals: 

,(a/3) _ [ J {a) (/3) (a/3) _ f J (a) ^x!^ 

For the ground HH state, we have 
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For both the LH-like and HH-like states, the z component of O is zero. It is important to note 
that f2 is proportional to the product / dz xd'^x' / dz"^ I dz xdx' / dz for the terms proportional 
to k and to the product / dzxx' I dzxdx'/dz for those proportional to k^, where x ^i-nd x' 
are the envelope functions obtained as a solution of Schrodinger equation with the potential 
U{z) and appropriate effective masses. For a symmetric U{z), the products of this kind are 
zero and the spin splitting vanishes in accordance with the general symmetry requirements. 

Based on these results for hole spectrum, we can calculate the related spin relaxation 
rate. Basically, two mechanisms of relaxation must be addressed. The first is the D'yakonov- 
Perel' mechanism.i^ It is related to spin precession with a frequency fl{kx, ky), which changes 
randomly due to rapid electron transitions in the momentum space. For this mechanism, 
the spin relaxation times for the x, y, and z spin components obey the following relation 
Tx = Ty = 2Tz = Top. For non-degenerate carriers, it is given as 
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In this equation, /o is the Boltzmann distribution function and r is an average characteristic 
time of electron scattering;— in the following, we will assume it to be equal to the electron 
momentum relaxation time. Another dominant source of spin relaxation is the Elliott- Yafet 



mechanism. It is caused by hybridization of the hole states at the finite wavevectors. Using 
the results for hole wavefunctions, an order-of-magnitude estimate in a strongly asymmetric 
QW is obtained for the ratio of Elhott-Yafet relaxation time Tgy to Top: 
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where Eq is the characteristic energy of hole confinement, typically of the order of few tens 
of meV. The expression suggests that the D'yakonov-Perel' mechanism provides the major 
contribution to spin relaxation in the case of well-defined quantized energy levels. Of course, 
this is not true for weakly asymmetric QWs, where the D'yakonov-Perel' mechanism can 
become less important. 

III. RESULTS AND DISUCSSION 

As a specific example, we consider SiGe inversion layers. The first case corresponds to a 
thin Ge layer grown on SiGe, and the second a Si layer on SiGe; an insulator is placed on 
top in each case. Such structures have been studied extensively due to their perspective ap- 
plications in p-type enhanced-mobility metal-oxide-semiconductor field-effect transistors.— 
Hole confinement in these devices is achieved by a strong perpendicular electric field and the 
insulator as schematically illustrated in Fig. 1. The barrier at the left-hand side is assumed 
to be infinitely high. The deformation energy in this case is = —2b i^cyijcw + 1) 5, where 
Cii and Ci2 are the elastic constants of the QW material (i.e.. Si or Ge), h is the deformation 
potential constant, and 6 is the relative mismatch in the lattice parameters. For the case 
of a relaxed SiGe buffer, the ground state in the Si QW is LH-like with E'^ > 0, while the 
Ge QW with E'^ < has a HH-like ground state. Note that the band discontinuity at the 
Si/SiGe or Ge/SiGe interface is not a well-defined quantity since not only the band gap but 
also the spin-orbit energy A are different across the interface. In our calculations, we assume 
that the confining electric field is strong enough to shift the holes away from the interface. 

In Figs. 3 and 4, we show the obtained spin splitting and spin relaxation rate l/T^p 
as a function of the confining electric field at T = 300 K. Figure 3 presents the case of a 
Si/Sio.yGeo.s structure, while Fig. 4 is for Ge/Sio.sGeo.y. The momentum relaxation time r 
of 5 X 10~^^ s is used, which corresponds roughly to the reported values of hole mobility . 



The material parameters are taken from Ref. UJ. The spin splitting is calculated for the 
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absolute value of the wavevector corresponding to the in-plane kinetic energy of ksT along 
the [100] (dotted line) or [110] (dashed line) direction, respectively. For the case of a Si QW, 
the main contribution to Q at room temperature is due to the terms proportional to k^. This 
is because the fc-linear term is proportional to B, which is small in the case of Si. On the 
other hand, the main contribution to Q in the Ge QW is from the HH— LH coupling since 
the spin-orbital gap A in Ge is large. Note that the mixing between valence and conduction 
band states can noticably influence spin splitting in the latter case (i.e., the Ge QW).— For 
a Si QW, this is less likely to be important since the direct band gap is 3.5 eV, much larger 
than that for Ge (0.9 eV).^^ 

As can be seen for the figures, the valence band of a Ge QW is characterized by spin 
splitting comparable to that for electrons in III-V QWs. Simultaneously, the hole spin 
relaxation time is quite short. In contrast, spin splitting and the relaxation rate for holes 
in a strained Si QW are much smaller. The reason for this is a weak spin-orbit interaction 
in Si, which is characteristic for light elements. In particular, the contribution to Q is 
proportional to A, which is as little as 44 meV for Si. Hence, it can be suggested that 
by applying a moderate longitudinal electric field, hole spin in a strained Si QW can be 
transferred over a distance of a micron or so without the loss of coherence. On the other 
hand, spin relaxation can be tuned effectively by modulating the confining electric field. 
It is important to note once more that in unstrained structures, say conventional p-type 
Si/Si02 inversion channels, spin relaxation is expected to be much faster due to relatively 
small subband separation and, as a result, strong hybridization of the hole states at finite 
wavevectors. 



IV. CONCLUSION 



We predict that spin splitting of 2D holes and the related spin relaxation time in strained 
SiGe QWs can be effectively controlled by the structure composition and the degree of QW 
asymmetry. For field-effect transistor-like structures,— spin splitting depends strongly on 
the value of the perpendicular electric field. In particular, the hole spin relaxation time in 
a strained Si QW channel can be in the hundred picosecond range at room temperature, 
making it suitable for spintronic applications. 
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FIG. 1: Model valence band profile used for spin splitting calculations. Holes are assumed to be 
confined in a triangular well with an infinite barrier on the left-hand side. The band discontinuity 
at the SiGe interface depends on the structure type (Si/SiGe or Ge/SiGe) and was disregarded. 
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FIG. 2: Spin splitting and D'yakonov-Perel' spin relaxation rate for a strained Si QW grown on 
the Sio.yGeo.s buffer. Spin splitting is given for the hole kinetic energy of ksT along the [100] 
(dotted line) and [110] (dashed line) directions. 
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FIG. 3: Spin splitting and D'yakonov-Perel' spin relaxation rate for a strained Ge QW grown 
on the Sio.3Geo.7 buffer. Spin splitting is given for the hole kinetic energy of fc^T along the [100] 
(dotted line) and [110] (dashed line) directions. 
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